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Abstract 

We derive the linear Langevin equation that describes the behavior of the fluctua- 
tions of the order parameter of the chiral phase transition above the critical tempera- 
ture by applying the projection operator method to the Nambu-Jona-Lasinio model 
at finite temperature and density. The Langevin equation relaxes exhibiting oscilla- 
tion, reveals thermalization and converges to the equilibrium state consistent with 
the mean- field approximation as time goes on. With the help of this Langevin equa- 
tion, we further investigate the relaxation of the critical fluctuations. The relaxation 
time of the critical fluctuations increases at speed as the temperature approaches 
toward the critical temperature because of the critical slowing down. The critical 
slowing down is enhanced as the chemical potential increases because of the Pauli 
blocking. Furthermore, we find another enhancement of the critical slowing down 
around the tricritical point. 

Key words: Chiral phase transition; Langevin equation; Critical slowing down; 

Mode coupling theory 

PACS:05.10.Gg,11.10.Wx,11.30.Rd 



1 Introduction 



It is widely believed that hadronic matters undergo a phase transition at high 
temperature, and hence become the quark-gluon plasma (QGP), where quarks 
and gluons are deconfined and chiral symmetry broken in the hadron phase is 
restored. Relativistic heavy-ion collisions provide significant opportunities to 
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explore the QGP. Heavy-ion collisions are essentially nonequilibrium processes 
and hence the description of the time evolution is necessary to understand 
the phenomena in a comprehensive way. Then, the hydrodynamic model is a 
powerful tool, and fairly correctly describes the experimental data of heavy- 
ion collisions, for example, transverse momentum dependence of the elliptic 
flow coefficient V2 [1,2,3]. 

Recently, our attention has been focused on dynamics near the critical points 
of the chiral phase transition [4,5,6,7,8,9,10,11,12,13,14,15,16,17,18]. Unfortu- 
nately, the hydrodynamic model (in particular, the so-called ideal fluid) will 
be inadequate to describe nonequilibrium phenomena near the critical point. 
Because it is valid only in the case where we are interested in gross variables 
associated with macroscopic time and length scales (entropy density, energy 
density, etc) and they are widely separated from other variables associated 
with microscopic scales(degrees of freedom of each quark and gluon). To un- 
derstand it, we should remember the dynamical hierarchy of classical dilute gas 
that is shown in TABLE 1 [19,20,21]. The classical gas consists of a lot of classi- 
cal particles and the dynamics is described by the Liouville equation. Then, the 
typical length and time scales are given by ro and ro/vth, respectively, where 
To is the interaction length and Vfh is the thermal averaged velocity. However, 
when we are interested in slow motions associated with the scales of the mean 
free path L and L/vth, the Liouville equation is reduced to the kinetic equation 
like the Boltzmann equation. Furthermore, when we observe quantities asso- 
ciated with further slower scales of the wave length and frequency of sound 
Lmacro ^nd T^acro, the Boltzmann equation is reduced to the hydrodynamic 
equation. Thus, in order to apply the hydrodynamic model, the microscopic 
scale must be smaller than the macroscopic scale, 6 — f'o/ Lmacro 1- a 
matter of fact, we assume the local equilibrium in the hydrodynamic equation 
and the deviation from the local equilibrium can be characterised by e. For 
example, the number density {n{x.,t))ini is approximated as 

{n(x,t))ini = {n(x))iocai + 0(e), 

where ( )ini means an expectation value by an initial state and ( )iocai denotes 
that by a local equilibrium state. When we completely ignore the deviation 
from the local equilibrium, we obtain the ideal fluid. When we take into ac- 
count the deviation from the local equilibrium, we obtain the dissipation term 
and the Euler equation is replaced by the Navier-Stokes equation [22]. 

Now, we return to the discussion of the critical dynamics. In the following, 
we assume that there exists a similar dynamical hierarchy even in quantum 
field theory because of the success of the hydrodynamic model in relativistic 
heavy ion collisions, although we do not have clear evidence so far [23]. On 
the other hand, to apply the hydrodynamic model for describing the critical 
dynamics, there must exist the clear separation between the microscopic scale 



2 



Region Equation length scale time scale 



microscopic 


Liouville 


^0 


ro/vth 


mesoscopic 


Boltzmann 


L 


L/vth 


macroscopic 


Hydrodynamic 


Lfnacro 


T 

-'-macro 



Table 1 

The dynamical hierarchy of classical dilute gas. vq is the interaction length, vth is 
the averaged thermal velocity and L is the mean fee path. L^acro f'lid T^acro mean 
macroscopic scales like wave length and frequency of sound. 

and the macroscopic scale, as is discussed above. However, near the critical 
points, the microscopic correlation length indefinitely increases and hence the 

expansion parameter e is not small any more. Therefore, we cannot apply the 
hydrodynamic model to describe the critical dynamics. On the other hand, it is 
probable that we still do not need the full details of the Heisenberg equation of 
motion because of the universality of phase transition. Then, from an analogy 
to gas dynamics, the critical dynamics can be regarded as a sort of mesoscopic 
scale dynamics and will be described by a kinetic equation. As a matter of 
fact, it is known that the critical dynamics in condensed matter physics is 
described by the Langevin equation. This theory is called the mode coupling 
theory [19,20,21,24,25]. 

One of the aims of this paper is to derive a linear Langevin equation for de- 
scribing the dynamics of the critical fluctuations above the critical tempera- 
ture of the chiral phase transition. There are two typical methods to derive the 
Langevin equation from a microscopic point of view [26] . One is the variational 
method that is applied to a sort of an effective action [27,28,29,30,31,32,33,34]. 
Then, the noise is introduced as an auxiliary field. However, as is discussed in 
[35] , there is arbitrariness for the introduction of the noise in this approach. 

In this paper, we use the projection operator method. This method was pro- 
posed by Nakajima [36], Zwanzig [37] and Mori [38], and developed from var- 
ious points of view [35,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54]. In the 
following, we will discuss based on the formulation evolved by Shibata, Hashit- 
sume, Uchiyama and the authors [44] . The derivation of the Langevin equation 
at vanishing chemical potential was discussed in [46] . In this paper, we extend 
the result to derive the Langevin equation at finite temperature and density. 
We show that our Langevin equation reveals thermalization and converges to 
the equilibrium state consistent with the mean-field approximation as time 
goes on. 

The other aim is to investigate the critical slowing down in the chiral phase 
transition [55] . As is well-known, the long wave length component of the order 
parameter shows extraordinary large fiuctuations near the critical point and 
the increase of the relaxation time. This slowing down of the relaxation of 



3 



the critical fluctuations is known as the critical slowing down. We define the 
relaxation time of the critical fiuctuations and investigate the temperature and 
chemical potential dependences. Then, we find the enhancement of the critical 
slowing down in the low temperature and large chemical potential region and 
around the tricritical point. 

This paper is organized as follows: In Sec. 2, we review the projection operator 
method following [44]. In Sec. 3, we introduce the Mori projection operator. In 
Sec. 4, we apply the projection operator method to the Nambu-Jona-Lasinio 
(NJL) model [56,57]. The equation obtained in Sec. 4 is not yet the Langevin 
equation. To derive the Langevin equation, we employ the renormalization of 
the memory function in Sec. 5. In Sec. 6, we calculate the power spectrum and 
discuss that our Langevin equation is consistent with the well-known results 
in equilibrium: the critical temperature, thermalization and the soft mode. In 
Sec. 7, we discuss the critical slowing down near the critical points by using 
the Langevin equation, and show that there exists the enhancement of the 
critical slowing down in low temperature and large chemical potential region 
and around the tricritical point. Concluding remarks are given in Sec. 8. 



2 Projection operator method 



In this section, we review the projection operator method. There are two dif- 
ferent approaches to derive coarse-grained equations in the projection operator 
method. When we apply this method in the Heisenberg picture, we can derive 
the Langevin equation. On the other hand, the master equation is derived 
in the Schrodinger picture [42]. In this paper, we discuss in the Heisenberg 
picture, following [44]. In the quantum field theory, the time evolution of op- 
erators is governed by the Heisenberg equation of motion. 



-0{t) = z[H,0{t)] (1) 
= iLO{t) (2) 
^O(i) = e^^(*-*°)0(to), (3) 

where L is the Liouville operator and is an initial time at which we prepare 
an initial state. The Heisenberg equation contains the information not only 
of gross variables but also of microscopic variables. The latter is irrelevant 
information and we carry out coarse-grainings by introducing a projection 
operator P. The projection operator P and its complementary operator Q — 
1 — P have the following general properties: 



PQ^QP^O. 



(5) 



From Eq. (3), one can see that the time dependence of operators is determined 
by e^H't--to)_ ^phis operator obeys the differential equation, 



at 

Prom this equation, we can derive the following two equations: 



±^iL{t-to)p ^ e'^it-t^) PiLP + e'^^*-^°'^QiLP, (7) 
dt 

^^iL{t-to)Q ^ e'^{t-to)piLQ + e'^^^-^°^QiLQ. (8) 
at 



Equation (8) can be solved for e*^^* *°^Q, 



iLit-to)Q ^ Q^iLQit-to) ^ J dre'^^^-^°^ PiLQe'^^^*-^\ (9) 



to 



Substituting Eq. (9) into Eq. (6) and operating 0{to) from the right, we obtain 
the time-convolution (TC) equation. 



d 



^^0{t)^e'^^'-'°^PiLO{to) + J dTe'^'^'-^^PiLQe'^'^^^-'°^iLO{to) 

to 

+Qei^Q(^-^o)iLO(to). (10) 



The first term on the r.h.s. of the equation is called the streaming term and 
corresponds to a collective oscillation such as plasma wave, spin wave and so 
on. The second term is the memory term that causes dissipation. The third 
term is the noise term. We can show that the memory term can be expressed 
as the time correlation of the noise. This relation is called the fluctuation- 
dissipation theorem of second kind (2nd F-D theorem). See Appendix A for 
details. The elimination of microscopic variables using the projection operator 
gives rise to the dissipation term and the noise term in the coarse-grained 
macroscopic equation. 

The TC equation is still equivalent to the Heisenberg equation and we cannot 
solve it exactly in general. In order to approximate the memory term, we 
reexpress it as 
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to 

where 



z 



t tl tn-l 



V{t, to) = l + Y, i"' Jdtijdt2--- J dtn 

"■=^ to to to 

xQLf{tn - to)QLf{tn-i -to)--- QLfih - to), 



(12) 



and 



Lf(t - to) = e^Q-^o«(*-*o)L/e-^^-^°^(*-*°). (13) 

Here, we introduced Lq and Lj, that are Liouville operators of the nonpertur- 
bative Hamiltonian Ho and the interaction Hamiltonian Hj, respectively, 



LoO^[Ho,0], LiO=[Hi,0]. (14) 
When we expand Vlt, to) up to first order in terms of Lj, we have 



d 



0{to) = e'^^'-'^^ PiLOito) + j dTe'^'^^-^^PiLQe"^^°^'^^-'°\LO{to) + ^{t). 

(15) 



dt 

to 



The noise term ^(t) should be determined so as to satisfy the 2nd F-D theo- 
rem, as we will see later. Equation (15) is the starting point in the following 
calculation. 



3 Introduction of the Mori projection 



In the derivation of the TC equation, we have not fixed a projection operator. 
If we can solve the TC equation exactly, the final result does not depend on the 
projection operator because the TC equation is equivalent to the Heisenberg 
equation. However, normally, it is impossible to solve the TC equation exactly, 
and we must start from the approximated equation (15). Then, the choice of 
the projection operator is important. There are several possible projection 
operators that extract slowly varying parts from an operator. In this study. 
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we adopt the Mori projection operator. The Mori projection operator projects 
any operators onto the space spanned by gross variables. There are three 
candidates for gross variables: (i) order parameters, (ii) density variables of 
conserved quantities and (iii) their products. When we can find out a complete 
set of gross variables, the macroscopic time evolution will be approximated by 
the superposition of the gross variables, as is shown in Fig. 1, schematically. 
For this purpose, we introduce the Mori projection operator defined by 



PO = ^qA, (16) 

i 

where Ai is the complete set of gross variables and the coefficient q is given 

by 

Q = ^(0,A,)-(A,A)T/. (17) 

j 

The inner product is the Kubo's canonical correlation, 

"rdX 

(X,r) = y -IV[p e^^Xe-^^y], (18) 
^ 

where p — e~^^/Tr[e~'^^] with the temperature (5"^. The inverse of the canon- 
ical correlation is defined by 

Y.i^,A)-^.{A,,A,)^5,,,. (19) 
j 

The physical meaning of the Mori projection operator is discussed in Appendix 
B in detail. 

4 Linear Langevin equation in the Nambu-Jona-Lasinio model 



We apply the projection operator method to the two-flavors and three-colors 
Nambu-Jona-Lasinio (NJL) model in the chiral limit and derive the linear 
Langevin equation for the critical fluctuations in the chiral phase transition at 
flnite temperature and density. The behaviors of the critical fluctuations are 
different between above and below the critical temperature Tc- In this study, 
we limit our discussion to the case oiT > T^.. 



The NJL Hamiltonian at flnite chemical potential is [56,57] 



7 



o 



5A2 



5Ai 



PO 



Fig. 1. The schematic figure of the Mori projection operator. An arbitrary operator 
O is projected onto the space spanned by gross variables 5Ai. The variable 5a means 
a microscopic variable. 



H = Hq + Hi, 

Ho = I ci3xg(x)H7- V-/X7°)g(x), 



(20) 
(21) 

(22) 



where Hq and Hj are the nonperturbative Hamiltonian and the interaction 
Hamiltonian, respectively. Here, we introduce the Pauli matrices r* {i = 
1, 2, 3). We set the coupling g = 5.01 GeV^^ and the three dimensional cutoff 
A = 650 MeV so as to reproduce the pion decay constant = 93 MeV and 
the chiral condensate (qq) = (—250 MeV)^ in the chiral limit. In particular, 
we are interested in the scalar channel. Thus, we ignore the pseudoscalar part 
of the interaction Hamiltonian in the following calculation. 

The quark fields q{x, t) and g(x, t) are expanded as 



g(x, t) = -= J:K,Mp, 5)e^P-e-^^^* + d^,MP, s)e-'^-e''^^% 



(23) 



(24) 



where Ep = |p|. The normalization conditions of the wave functions u{p,s) 
and v{p, s) are 



u^{p, s)u{p, s') = v\p, s)v{p, s') = Ss,s', 
v{p, s)u{p, s) = v^{p, s)u{-p, -s') = 0. 



(25) 
(26) 
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The corresponding commutation relations are given by 



[g(x,t),g(y,t)]+ = 7V3nx-y), 

[^(x, t), g(y, t)]+ = [g(x, t), g(y, t)]+ = 0, 



(27) 

(28) 
(29) 



where [ ]+ means the anticommutator. 

Now, we choose gross variables to define the Mori projection operator. As is 
discussed in Sec. 3, there are three possible candidates for gross variables: (i) 
order parameters (ii) density variables associated with conserved quantities 
and (iii) their products. For simplicity, we ignore (iii) in this study. Then, 
we have the order parameter qq, the number density q'j^q, the energy density 
T'^^ and the momentum density as gross variables from the conditions 
(i) and (ii). However, we assume that near the critical temperature of the 
phase transition, the order parameters are much slower than the other density 
variables because of the critical slowing down where the relaxation time of the 
critical fluctuations indeflnitely increases (At high density, the number density 
can be another gross variable as in the glass transition. We will discuss this 
point in Sec. 8.). Thus, we exclude (ii) from our gross variables. After all, the 
gross variable relevant to our calculation is the order parameter of the chiral 
phase transition. Then, the Mori projection operator is 



where (g(x)g(x))e5 is the expectation value in thermal equilibrium and van- 
ishes above the critical temperature [58]. The canonical correlation is 




(30) 



Here, we have used the notation. 



5a{x) = ?(x)5(x) - (g(x)g(x)) 



(31) 




(32) 



where 



{O)o 



^0 



(33) 
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with Zq = Tr[e~^^°]. If we can solve the equation exactly, the system should 
thermalize with the total Hamiltonian H. However, in the following calcu- 
lation, we make approximation where the solution of the derived Langevin 
equation relaxes toward an equilibrium state consistent with the mean-field 
approximation where quarks behave massless free particles above the critical 
temperature. Thus, we replace H with Hq in the definition of the canonical 
correlation (18). In this manner, we should choose the projection operator so 
as to abstract the process that we should describe. In this sense, the projec- 
tion operator method has more than the lowest order approximation of the 
perturbative expansion. 

Substituting this Mori projection operator into the TC equation (15), we can 
derive the Langevin equation of the critical fluctuations of the order parameter 
in the chiral phase transition. 



(Jjv j 



+e(x,i). (34) 

The noise (third) term is determined so as to satisfy the 2nd F-D theorem. 
And, as is shown in Appendix C, the streaming term vanishes. Thus, in the 
following, we discuss the calculation of the memory function. 

The calculation of the memory term is quite intricate because we must cal- 
culate the coarse-grained time evolution operator e^'^^oQt instead of the nor- 
mal time evolution operator e*^°* that appears in normal calculations without 
coarse-grainings. First, we expand the memory term with the help of the nor- 
mal time evolution operator: 



^ dTe'^^'-''^PiLQe'^^°'^^iL5a{x.) = J dre'^^'-^^ PiLe'^°^B{T,0)QiL5a{x), 



(35) 



where 



tn-l 



B{t, to)^l+ j dh--- j dtnL^ih - to)U{t2 -to)--- U{tn - h), 

to to 

(36) 



with 
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(37) 



Each term does not have the coarse-grained time evolution operator and is 
calculated exactly. Finally, we find the n-th order term given by 



J dre'^^^-^^PiLe'^^^i-i)'' J dti--- dtnU{ti) ■ ■ ■ L^{tn)QiL5o{^) 



t T 

= (-1)"/ dre'^^^-^^PiL J dh ■■■dtn 



d^Xt„ (x, x„) dXt^_^_t^ (x ) (iXt^_t2(x2,xi) 
^ 1^, dF, Mxi,r-tO, 

(38) 

where (5cro(x, t) — e'^^°*5a{x,0) and 



Xt(x,x") = (5c7o(x,t),5a(x')) • (5a(x'),5a(x"))-' 

= |ciVxi(x-x')Xo^(x'-x"). (39) 

In this expression, we dropped the sign of the integral for the repeated space 
variables for simple notation. In the following, we use this notation without 
notice. The concrete form of Xt(x) is given later. Thus, the memory term is 
given by 



t 

J dTe'^^^-^^PiLe'^^^B{T, Q)QiL6a{^) 



* j2 j oo * 

= j dri—Xri^, xi) • 5a(xi, t-T)+ e^^(*--)PiL,— 5<7o(x, r)) + YlJ M-^T 



X / d^xidii • • • d^x ^^^^n(x,x„) _ tiXt„_^-t„(x„,x„_i) _ _ _ dXt,-t^{^2,y^i) 
y " " dtn-i dti 

•(^X,_t,(xi, y) • 5a(y, t - r) + e^^(*-^)PiL,5(7o(xi, r - ii)). (40) 
Here, we used e'^(*~'^Mcr(x) = 5a{x, t — r). 

This expression is simplified by using the Laplace transform. First, we ignore 
the terms including the interaction iL^. The Laplace transform of the first 
term is 
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£{— X4(x,xi) ® (5a(xi,t)} = (x,xi) ■ (5(7^(xi,s), 



where C{ } means the Laplace transform and 



(41) 



C{A ®B}^C{j dTA{t - t)B{t)}. (42) 



Here, we introduced the following expressions, 

oo 

£{(7(x, t)} = j die-'V(x, t) = (7^(x, s), (43) 



£{-X,(x,Xi)} = / die-*-X,(x,Xi) = X,^(x,xi), (44) 



£{— X,(x,xi)} = / (iie-^*^X,(x,xi) =Xf(x,xi). (45) 



The Laplace transform of the second term is 



t 



(fXtJ^ xi) dXt-t^ (xi, y) 
dtl 

KJ 

= -X,^(x,xi) •Xf(xi,y) •5a^(y,s). (46) 



c{-jdt, — ^ 0My,t)} 



All other terms can be calculated in the same way. Thus, the sum of all the 
Laplace transform of the memory term is given by 



* ,2 oo * 

C{J ds—X,{^,^i)-6a{^^,t-s) + J2i-'^T I dr J dh---dtn 



Here, we used that 



X*(k) 
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-J d^(x - x')e-^^("-"')(5(7o(x, t), 5a{^', 0)) 



^ EpEp+k - p(p + k) ^i(E^-E^^^)t 

EpEp+]^{Ep — -Bp+k) 

, . . ^p^p+k + p(p + k) 

-C^p-C'p+kl.-C/p + -C/p+kJ 

X (e'(^p+^p+'')* + e"^(^p+-^p+'')*)}, (48) 
where = 3, Nf = 2, E^ = |k| and 



= , . (49) 

The Laplace transforms are defined by 



(k) = I dte-^'j^xm, (50) 



e(k) = / dte-^*— x*(k). (51) 



This term has the well-known structure and can be interpreted in the kine- 
matical way. The first term is the scattering process and the second term is 
pair creation and annihilation processes. 

Next, we calculate the terms including iL/. The uncalculating factor PiLi6ao{x, t) 
in Eq. (40) contains various contributions. Here, we are concerned with the 
contribution of the particle- antiparticle loop diagram shown in Fig. 2. (Then, 
our Langevin equation shows consistent behaviors with the equilibrium results 
in the mean- field approximation, as we will see later.) In this random phase 
approximation, we can calculate the interaction part by using Wick's theorem. 



(iL/5cro(x,r),(5a(z)) 

~ -ig I ciV2[ gg(y, -zA)g ^(x, r - iA) gg(y, -iA)g o^(x, r - iX) ql{y, -iX)f {z) q^{y, -iX)q '{z) 

- go'(x,r-iA)^ g(y, -iA) g(]'(x, r - iA)| ^(y, -iX) q^{y, -iX)q ^{z) gg(y, -iA)g ^^(z) ] 
= -2g/3 1 dV^Xr(x, y)xo(y, z). (52) 

In the second line of the r. h. s. of the equation, we approximate the function 
by calculating the contribution of the ring diagram shown in Fig. 2. Here, 
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Fig. 2. The ring diagram that contributes the calculation of the memory term. 

we introduced the contraction defined in Appendix D. Finally, the Laplace 
transform of the remaining part of the memory term is given by, 



n n=l n 





X dti---dtn e ^ 'PiLidao{xi,T -ti))} 

J dti dti 

" 

(53) 

Summarizing the above results, the memory term is given by 



J dTe'^^^-^^PiLQe'^^°^^iL5a{x) ^ - J dr J dVr(x - x', t)(5(7(x', i - r). 



(54) 

Here, the memory function r(x, t) is given by the inverse-Laplace transform 

of r^(x,s). 



(66) 



In short, the memory function that includes the coarse-grained time evolution 
operator e*'^'^"'?* is expressed by the combination of the normal correlation 
function Xt(k). This fact was pointed out in [59] for the first time. 

Finally, we discuss the noise term. In principle, the noise term also can be 

calculated by using the perturbative expansion. However, we do not need the 
concrete form of the noise in the following calculations. Thus, we determine 
the noise so as to reproduce the 2nd F-D theorem. Then, the equation is given 

by 
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Fig. 3. The time dependence of the memory function at ^ = 200 MeV. The sohd, 
dashed and dotted hnes are the memory function for temperatures T — = 10 
MeV, 30 MeV and 80 MeV. 



j^5a{k, t)^-j cirr(k, r)5a(k, t - r) + ^(k, t), 



(56) 



where the introduced noise ^ (k, t) has the first and second order correlations 
given by 



(e(k,t),Mk',0)) = (e(k,t))o = 0, 
(e(k, i), e(k', t')) = V5tlT% t - Oxo(k). 



(57) 
(58) 



The above equation results from the naive application of the projection oper- 
ator method. However, we cannot interpret this equation as a Langevin equa- 
tion. To derive the Langevin equation, we must employ the renormalization 
of the memory function, as we will discuss in Sec. 5. 



5 Renormalized Langevin equation 



We have applied the projection operator method to the NJL model and derived 
the equation (56). Unfortunately, we cannot regard this equation as a Langevin 
equation, because of the problem of the long time correlation of the memory 
function. In Fig. 3, the time dependence of the memory function r(0,t) is 
shown at /X = 200 MeV. The sohd, dashed and dotted hnes are the memory 
function for temperatures T—T^ = 10 MeV, 30 MeV and 80 MeV. The memory 
functions converge to finite values at late time. This behavior is unphysical 
because the memory function is given by the time correlation function of the 
noise from the 2nd F-D theorem. It is clear from the definition of the noise (10) 
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that because of the projection operator Q, ah gross variables are excluded from 
the time evolution of the noise and hence the time correlation must converge 
to zero rapidly. 

The long time correlation means that a variable associated with a long time 
scale is still included in the noise because of the incompleteness of the definition 
of the Mori projection operator [61]. As an example, we consider an exactly 
solvable model (See Appendix E for details.). The Hamiltonian is 

It is known that we must choose two gross variables p and x to derive the 
correct Langevin equation. Then, we have 



-p{t) = -Mulx{t) - I dsE{t - s)p{s) + fit), 



(61) 

where the memory function E{t) and the noise f{t) are given by 



E{t) = if it), /(O)) • {p,p)-' = ^ E ^ ^osu;,{t - t'), 

^ j=i 



N ■ r 

= E%[cosa;,i(x,(0) - ^x(O)) + ^^p,(0)], 

3=1 



(62) 



(63) 



respectively. On the other hand, when we choose only one gross variable p, 
the resulting equation has the long time correlation. 



d 
Jt 



p(t) ^-j dsE'it - s)p{s) + fit), (64) 



where 



s'(t) = (fit), rm ■ ip,p)-' = uji + s(t), (65) 

f\t) = -Mujlx{Q) + f(t). (66) 
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Model 


GVl 


GV2 


MFl 


MF2 


F 


Exact 


P 


X 


E'{t) 


m 




NJL 


6a 


? 


T{k,t) 


$(k,i) 


J)2(k, t) 



Table 2 

The correspondence between the exact model and the Langevin equation in the NJL 
model. The first two columns indicate the gross variables (GV) in each model. The 
next two columns denote the memory function with long time correlation (MFl) 
and that without long time correlation (MF2). The last column shows the frequency 
function (F). 



T-T. = 10 MeV 
T-T„ = 30 MeV 
T-T- = 80 MeV 



-0.02 ' 1 1 1 1 1 

1 2 3 4 5 

t[fml 

Fig. 4. The time dependence of the renormalized memory function at /x = 200 
MeV. The solid, dashed and dotted lines are the memory function for temperatures 
r - Tc = 10 MeV, 30 MeV and 80 MeV. 



As is discussed in Appendix E, if we observe the system with long time scale, 
converges to zero quickly and does not have long time correlation. How- 
ever, does not converge to zero because of the frequency lJq and has 
long time correlation. As just described, the behavior of the memory function 
strongly depends on the definition of the Mori projection operator. In order 
to obtain the memory function without the long time correlation, we must 
prepare the complete set of gross variables. 

The correspondence of the exact model and our equation is given in TABLE 2. 
The gross variable p in the exact model corresponds to 5a(k., t) in our Langevin 
equation. Therefore, in order to derive the Langevin equation without long 
time correlation, we must find out another variable corresponding to x in our 
model and redefine the Mori projection operator with the two gross variables. 
However, another gross variable in the NJL model is not obvious because the 
NJL Hamiltonian is composed of q and q, and does not include the composite 
operator a — qq explicitly. 

Therefore, in this paper, we derive the correct memory function in another way. 
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Fig. 5. The time dependence of the frequency function at fi = 200 MeV. The sohd, 
dashed and dotted hnes are the memory function for temperatures T — = 10 
MeV, 30 MeV and 80 MeV. 

First, we assume that there exists the exact correspondence of the exact model 
with our equation, and speculate the correct form of the memory function. By 
comparison Eq. (64) with Eq. (56), our memory function corresponds to H'(t), 
while the correct memory function should be given by the counterpart of 
S(i). In order to obtain the correct memory function in our model, we should 
separate the r(k, t) into two parts: the frequency part and the dissipation 
part. In quantum field theory, the frequency is given by the real part of the 
self energy in energy-momentum space, that is, the imaginary part of the 
memory function. Thus, we separate the imaginary part from the memory 
function to redefine the memory function without the long time tail. As a 
result, our memory function is separated into the following two parts. 



The memory function without long time correlation is given by the "renor- 
mahzed" memory function $(k, t), while ^l^{t) gives the frequency function. 

The time dependence of the renormalized memory function and the frequency 
function a,t pL — 200 MeV are given in Figs. 4 and 5, respectively. The sohd, 
dashed and dotted hnes are the memory function for temperatures T — T^ — 10 
MeV, 30 MeV and 80 MeV. The renormalized memory function relaxes rapidly 
and vanishes at late time, while the frequency function converges to a finite 
value depending on the temperature. In this manner, we can get rid of the 
long time correlation of the memory function by renormalization. 




(67) 
(68) 



(69) 
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Consequently, the correlation properties of the noise arc also changed. The 
correct noise is related not to the memory function but to the renormalized 
memory function by the 2nd F-D theorem. At last, we have obtained the 
renormalized linear Langevin equation. 



;5a{k, t)^-J dTnl{T)5a{k, t - r) - j dT$(k, T)(5(7(k, t - t) + ^(k, t), 



(70) 

where 

(e(k,t),Mk',0)) = (e(k,t)) = 0, (71) 
(e(k, t), et(k', t')) = l^41$(k, t - t')xo(k). (72) 

From the behavior of the renormalized memory function, we can see the typical 
time scale of the noise is given by about 1 fm, as is shown in Fig. 4. As is 
discussed above, the time scale of the renormalized memory function is nothing 
but that of the noise. Thus, the typical time scale of the microscopic variables 
we have coarse-grained is about 1 fm. In our Langevin dynamics approach, 
the time scale of gross variables must be larger than that of the coarse-grained 
variables. Thus, our Langevin equation is reliable when the relaxation time of 
Sa{k. t) is larger than 1 fm. This condition is always satisfied near the critical 
points because of the critical slowing down. 

This Langevin equation is the equation of operators, that is, the quantum 
Langevin equation. However, in the following discussions, we solve Eq. (70) as 
a semiclassical equation with a noise as a classical random field, because the 
concept of the noise as an operator is not clear. For this purpose, we introduce 
a classical noise that reproduces the correlation properties defined above, 

< e(k, t)(5a(k', 0) > = < e(k, t) >= 0, (73) 
«^(k,t)r(k',t') »=(e(k,t),enk',t')), (74) 

where ^ ^ means the average for the noise with a suitable stochastic 
weight. However, even after this replacement, the quantum effect is still in- 
cluded in fl'^{k,t), $(k, t) and Xo(k). In this sense, this is a semiclassical 
Langevin equation. 

The correlations of the noise obtained by applying renormalization is consis- 
tent with those in the generalized Langevin equation for the Brownian motion 
discussed in [60], where the correlation of the noise is given by the real part 
of the memory function rather than the memory function itself as is discussed 
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Fig. 6. The time evolution of <^ 6a{0,t) S> at /i = 200 MeV. The soUd, dashed and 
dotted hnes are the <C 6a{0, t) > for temperatures T - Tc = 10 MeV, 30 MeV, 80 
MeV. 

above. Otherwise, the Langevin equation does not show thermahzation at late 
time. 

The averaged time evolution of the critical fluctuations of the order parameter 
at vanishing momentum is shown in Fig. 6. One can see that the nonequi- 
librium fluctuations relax exhibiting oscillation and finally converge to zero. 
This indicates that the critical dynamics of the chiral phase transition may 
not be described by a simple diffusion-type equation like the time-dependent 
Ginzburg-Landau (TDGL) equation. The relaxation time of the critical fluc- 
tuations becomes larger as the temperature approaches toward the critical 
temperature because of the critical slowing down that is investigated carefully 
in Sec. 7. 



6 Thermalization and power spectrum 

We have derived the linear Langevin equation in the projection operator 
method, employing the lowest order approximation of the perturbative ex- 
pansion and the renormalization of the memory function. In this section, wc 
show that the Langevin equation converges to an equilibrium state consistent 
with the mean-field approximation as time goes on. 

The relaxation time of the critical fluctuations diverges at the critical point 
because of the critical slowing down. Thus, the critical temperature of the 
second order phase transition can be determined by the temperature where 
$(0, t), and hence r(0, t), vanishes. From Eq. (55), the vanishing point is given 
by the following condition. 
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l-2^/3Xo(O)|T=Te = 0. 



(75) 



This condition is nothing but the self-consistency condition of the chiral con- 
densate calculated in the mean- field approximation [62] . Therefore, the critical 
temperature of second order phase transition in the linear Langevin equation 
is completely same as the conventional mean-field result. Besides, this condi- 
tion can be regarded as a kind of the Thouless criterion [63]. On the other 
hand, it is not clear that the critical temperature of the first order phase tran- 
sition can be determined by the same condition. However, from the numerical 
result discussed below, one can see that the critical slowing down occurs even 
for the first order phase transition. 

Next, we show that the linear Langevin equation reveals thermalization. From 
Fig. 6, 6a(k., t) converges to zero at late time for arbitrary initial states. This 
is one evidence of thermalization. The other evidence is the correlation func- 
tion of Sa(k,t). From the renormalized Langevin equation, we can derive the 
equation of the correlation function. 



- < 5a{k, t)5a{-k, 0) > = - y drnKr) < 5a{k, t - T)da{-k, 0) > 



t 

- J dr^k, t) < 5(T{k, t - T)5a{-k, 0) > . 



(76) 

When, we prepare the thermal equilibrium state as an initial state, the equa- 
tion represents the correlation function in equilibrium. 



< da{k, 0)Sa{-k, 0) >= {S&{k, 0), Sa{-k, 0)) = xo(k). (77) 

The solution of the differential equation is defined in t > 0. Thus, we assume 
that the solution is symmetric att — O. Then, the correlation function is given 

by 



< Sa(k,t)6a{-k,0) > 

2X0 (k) 



oo 




/ 


— Re 




2n 


— oo 




oo 


duo 




/ 

— oo 


27r 1 - 



-iujt 



-iuj l^k(^) + *(k,a;)_ 
2$(k,cu)Xo(k) 



27r I -ia; + f]^(cu) + $(k,cu)|2 



-ibjt 



(78) 



where 
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Fig. 7. The temperature dependence of the power spectrum at k = for a fixed 
chemical potentials = 0. The solid, dashed and dotted lines are for T — Tg = 1 
MeV, 10 MeV and 30 MeV, respectively. 



OO 

nl{uj) = J dml{t)e"^^ = ilmr^(k, -iu + e), 



oo 

$(k, u;)^J dt^k, t)e''^^ = Rer^(k, -iu + e). 



(79) 
(80) 



It should be noted that ilk(<^) $(k, a;) are the odd and even functions of 
cu, respectively. 

The same result can be obtained from the Wiener- Khinchin theorem, where 
the correlation function in equilibrium is given by 

lim < 6a{x + x', t + f)6a(x', t') > 

t'— >oo 



/(k,a;)e-*'"*e 



iujt 



(81) 



Here, the power spectrum /(k, a;) is defined by 



I(k,uj) = lim ^ \Scr{k,uj)\'^ >, 



(82) 



where 



1 

Sa{k,uj) = lim J dt6a{k,t)e'''\ 



(83) 
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Here, T is the period of the time evolution of (5(j(k, t). In this derivation, we do 
not use the information of the initial state. However, the result is completely 
same as Eq. (78). By this means, we obtain the same expression of the corre- 
lation function in different two ways under the assumption of thermalization. 
Furthermore, if we ignore the oscillation term, the above result reproduces the 
same result as the power spectrum in the Brownian motion [60]. In this sense, 
our Langevin equation reveals thermalization. 

And this is the reason that we adapted the renormalization procedure dis- 
cussed in the previous section. To get rid of the long time correlation of the 
memory function, it is possible to apply another condition of renormalization; 



However, it is clear that we cannot show thermalization in this renormalization 
condition. Therefore, in order to assure thermalization, we should employ the 
renormalization condition discussed in the previous section. 

The temperature dependence of the power spectrum is shown in Fig. 7. We 
can see that the peak moving toward origin becomes prominent as the temper- 
ature is lowered toward T^. The power spectrum characterizes the space-time 
correlation in the energy-momentum space and hence can be interpreted as 
the spectral function in thermal field theory. Then the peak with narrow width 
reveals the existence of a collective mode whose energy tends to vanish as the 
temperature approaches toward T^. Such a mode is called a soft mode. The 
soft mode appears when the system becomes unstable for external perturba- 
tions [64,65]. The temperature dependence of the power spectrum is consistent 
with the the previous result, where the spectral function is calculated in the 
linear response theory [64]. 



7 Enhancement of critical slowing down 

As is well-known, the long wave length component of the order parameter 
shows extraordinary large fluctuations near the critical point because of the 
increase of the relaxation time. This slowing down of the fluctuations is called 
the critical slowing down [55]. In this section, we discuss the behavior of the 
critical slowing down in T-fi plane. 

For this purpose, flrst, we calculate the phase structure of the NJL model. As is 
discussed, the renormalized Langevin equation (70) has the same critical tem- 
perature as that in the mean-field approximation. Thus, we should calculate 
the thermodynamic potential in the mean-field approximation in order to fix 




(84) 
(85) 



$(k,i) 



r(k,i) -r(k,oo). 
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Fig. 8. The phase diagram in T — plane in the mean-field approximation. The 
dashed and solid lines denote the critical line of the second and the first order phase 
transition, respectively. The tricritical point is located at (T, /x) = (70, 280) MeV. 
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Fig. 9. The thermodynamic potential per unit volume on the critical line of the 
chiral phase transition. The solid, dashed, dotted and dotted dashed lines represent 
the thermodynamic potential at /x = 0, 150, 280, 300 MeV, respectively. 

the phase diagram, that is shown in Fig. 8. The dashed line represents the crit- 
ical line of the second order phase transition and starts from (T, /x) = (185, 0) 
MeV. The order of the phase transition changes from second to first at the 
tricritical point located at (T, fi) — (70, 280) MeV. The solid hne expresses the 
critical line of the first order phase transition that ends at (T, /i) — (0, 317) 
MeV. 

The corresponding thermodynamic potential per unit volume on the critical 
hne is shown in Fig. 9. The solid, dashed, dotted and dotted dashed fines rep- 
resent the thermodynamic potential at = 0, 150, 280, 300 MeV, respectively. 
The bottom of the thermodynamic potential becomes broader as the chemical 
potential approaches toward the tricritical point. This reflects the existence 
of the large fluctuations around the tricritical point. At — 300 MeV, the 
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Fig. 10. The temperature dependenee of the power spectrum at k = for a fixed 
chemical potentials jj. = 280 MeV. The solid, dashed and dotted fines are for 
T-Tc = l MeV, 10 MeV and 30 MeV, respectively. 

order of the phase transition is first and hence the thermodynamic potential 
has two local minima, although it is invisible. 

In order to clarify tlic equilibrium property above the critical temperatures, 
we calculated the temperature dependences of the power spectrum at // = 280 
MeV in Fig. 10. Compared to Fig. 7, one can see that above the second 
order phase transition, the temperature dependence of the power spectrum 
is almost same, although the peaks at /i = 280 McV approach to the origin 
faster. However, interestingly enough, the behavior of the memory function 
shows that the relaxation is drastically changed around the tricritical point as 
we will see soon later. 

We are interested in the behavior near the critical points and the typical time 
scale of the critical fluctuations increases because of the critical slowing down. 
Then, we can ignore the microscopic time-dependence included in the memory 
function and the frequency function. As is shown in Fig. 4, the rcnormalized 
memory function converges to zero with a short time scale about 1 fm, that 
is negligible by comparison with the time scale of the critical fluctuations 
in the vicinity of the critical temperature. Then, we can approximate the 
renormalized memory function as 




(86) 



where 



oo 




(87) 
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On the other hand, the frequency function converges to a finite value rapidly 
and we approximate it as 



Q^(k, t)^nl= lim Q^(k, t). (88) 



t— >oo 



These approximations are sometimes called the Markov limit. Then, the Marko- 
vian Langevin equation is given by 



f^Sa{k, t) = -nl J drSaik, r) - 7k5a(k, t) + ^(k, t), 



(89) 

where 

< C(k, t)5a(k', 0) > = < e(k, t) >= 0, (90) 
« e(k, m^', t') » = 2y7kXo(k)4%<^(^ - t'). 

(91) 

In the Markov limit, the colored noise is reduced to the Gaussian white noise. 
There are two possibility of the solution of the Langevin equation, depending 



on the value of = y^^k ~7k/4- The critical fluctuations relax exhibiting 
oscillation, if A\ > 0. However, if < 0, this is the case of the overdamping 
and the solution does not show oscillation. In our case, the condition Al_> 
is always satisfied. Thus, the solution of this Langevin equation is 



•C Sa(k,t) > 



da{k,Q)e 



-lkt/2 



'^^ sin A]st + cos Ai^t 



2A 



(92) 



Prom this analytic expression, one can see that the typical relaxation time 
of the long wave length component of the critical fiuctuations is given by 

Trt = 2/70. 

First, we discuss the behavior of the frequency that is shown in Fig. 11 
at a vanishing momentum. Above the critical line of the second order phase 
transition, the frequency decreases at the almost same speed as the temper- 
ature approaches to the critical temperature for different fi. On the other 
hand, above the critical line of the first order phase transition, the frequency 
increases with /x along the critical line from the tri critical point. This behav- 
ior is consistent with the behavior of the thermodynamic potential. Above 
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Fig. 11. Temperature and chemical potential dependences of the frequency Aq. 
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Fig. 12. Temperature and chemical potential dependences of the relaxation time at 
higher temperature than 50 MeV. 



the second order phase transition, the thermodynamic potential has only one 
minimum and the curvature of the potential is negligibly small as is shown in 
Fig. 9. However, the thermodynamic potential has two minima above the first 
order phase transition and the curvature of the potential can be finite. 

The behavior of the relaxation is extremely interesting in comparison with that 
of the frequency. The temperature and chemical potential dependences of the 
relaxation time is shown in Fig. 12. The relaxation time at lower temperature 
than 50 MeV is not shown because it diverges at speed. Thus, we plotted 
the damping coefficient 70 in low temperature and large chemical potential 
region in Fig. 13 instead of the relaxation time. At a fixed /x, the relaxation 
time increases as the temperature is lowered toward the critical temperature 
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and diverges at the critical point because of the critical slowing down. In this 
figure, we started to plot from the 1 MeV higher temperature than the critical 
temperature. When we increase /x along the critical line, we encounter the 
enhancement of the critical slowing down; one is in the low temperature and 
large chemical potential region and the other is around the tricritical point. 
The former can be explained by the Pauli blocking due to the Fermi surface. 
As a matter of fact, the damping coefficient vanishes quickly there as is shown 
in Fig. 13. 

More interestingly, the critical slowing down can be enhanced around the tri- 
critical point. This would be related to the broadness of the bottom of the 
thermodynamic potential. As is shown in Fig. 9, when we increase the chemi- 
cal potential, the bottom of the thermodynamic potential increases gradually 
and achieves maximum at the tricritical point. At larger chemical potential 
than that of the tricritical point, the thermodynamic potential has two mini- 
mum but the potential barrier between the two minima is still shallow. Thus, 
the large fluctuations are still possible to survive. In short, there exist ex- 
traordinary large fluctuations around the tricritical point and this causes the 
enhancement. Because of the large relaxation time of the critical fluctuations, 
thermalization is decelerated, in particular, in the low temperature and large 
chemical potential region and around the tricritical point. 

It should be noted that the enhancement is not obvious in the behavior of 
the power spectrum. This is because the power spectrum is given by $(k, cu) 
and fl^{u}), and the former is smaller than the latter. Thus, the behavior of 
the power spectrum is dominated by ^l^{uj) and the information of <l>(k, cj) is 
smeared. However, as is indirectly shown in Fig. 11, the frequency does 
not reveal the anomalous behavior around the tricritical point, clearly. Thus, 
we cannot see the remarkable difference between Figs. 7 and 10. 



8 Concluding remcirks 

We have derived the linear Langevin equation that describes the dynamics 
of the chiral phase transition at finite temperature and density. The simple 
application of the projection operator method to the Nambu-Jona-Lasinio 
(NJL) model caused the long time correlation in the memory function. To 
avoid this difficulty, we introduced the renormalization of the memory func- 
tion. The resulting Langevin equation reveals the critical slowing down at 
the same temperature as that estimated in the mean-field approximation and 
shows thermalization. 

The order parameter relaxes exhibiting oscillation and this is different from 
the behavior expected from a simple diffusion-type equation like the time- 
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Fig. 13. The temperature and chemical potential dependences of the damping coef- 
ficient 7o in low temperature and large chemical potential region. 

dependent Ginzburg-Landau (TDGL) equation . This means that the assump- 
tion employed to derive the TDGL equation may lose its validity in the chiral 
phase transition. In the derivation of the TDGL equation, one assumes that 
the time evolution of the order parameter is induced by a sort of the ther- 
modynamic force that is given by 5F/SM, where F is the Ginzburg-Landau 
free energy and M is the order parameter. However, in general, a number of 
thermodynamic forces are possible to induce such an irreversible process be- 
cause of the cross effect ( for example, the Dufour effect and the Serot effect 
in thermal conduction). Then, the usual TDGL equation can be changed from 
a simple diffusion-type equation. 

The power spectrum was also calculated in different two ways. One is the appli- 
cation of the Wiener- Khinchin theorem and the other is to solve the evolution 
equation of the correlation function under the assumption of initial thermal 
equilibrium. The two power spectrums have completely same form and we can 
conclude that the Langevin equation shows thermalization. The energy and 
temperature dependences of the power spectrum is consistent with that cal- 
culated in the linear response theory where the random phase approximation 
is employed [64] and there exists the soft mode. As a result, we can conclude 
that our Langevin equation fulfils the requirements near Tp. 

In this calculation, we approximate the memory function calculating the ring 
diagram contribution because we want to derive the equation that converges 
to the equilibrium state consistent with the mean-field approximation. How- 
ever, it is possible to derive the equation including higher order contributions. 
Then, the merit of our formulation is to be possible to derive nonperturba- 
tive results by approximating the memory function perturbatively, and the 
higher order term is already known by the systematic expansion oiV{t,to). 



29 



In the latter part of this paper, we discussed the temperature and chemical 
potential dependences of the critical fluctuations in the chiral phase transition. 
At a fixed /i, the relaxation time of the critical fluctuations increases because of 
the critical slowing down. On the other hand, increasing the chemical potential 
along the critical line, we found the enhancement of the critical slowing down. 
One is seen in low temperature and large chemical potential region. This is 
because the decay of the critical fluctuations is suppressed owing to the Pauli 
blocking. The other enhancement is shown around the tricritical point. This 
would be related to the broadness of the thermodynamic potential around 
the tricritical point and means that there exist extremely large fluctuations. 
The enhancement of the critical slowing down may affect the formation of the 
disoriented chiral condensate (DCC). Around the tricritical point, the large 
fluctuations can have a long relaxation time and it may be of advantage to 
form the large domain of DCC. 

As is well-known, the mean-fleld approximation cannot describe the static 
and dynamical critical exponents correctly in various cases. Thus, the validity 
of our results might be suspected. However, it should be noted that we ap- 
proximate the memory function by the ring diagram contribution and hence 
the relaxation dynamics itself includes the fluctuation effects, although the 
equilibrium state is consistent with the mean-fleld one. 

As another problem, we used the NJL model as a low-energy effective the- 
ory of QCD, where the endpoint is located at lower temperature and higher 
chemical potential than those estimated in other approaches, for example. 
Lattice QCD calculations [66,67,68,69,70,71,72]. However, we can expect that 
our calculation is still available to describe the critical dynamics qualitatively, 
because the enhancement of the critical slowing down can be explained from 
the Pauli blocking and the behavior of the thermodynamic potential. Thus, 
it is reasonable to expect that the behavior still survive even if we take into 
account higher order corrections or calculate in more realistic models. 

In this study, we ignored the finite current quark mass. When we take into 
account it, the phase transition is changed from the second order to crossover 
and the tricritical point is replaced by the endpoint. Then, the critical slowing 
down is smeared and hence the behavior around the endpoint will be changed. 

We have discussed the derivation and the behavior of the Langevin equation 
ignoring nonlinear effects. However, from the Ginzburg criterion, the nonlinear 
fluctuations play important roles near the critical points. As a matter of fact, 
the mode coupling theory makes it clear that to calculate the dynamical critical 
exponent correctly, we must introduce nonlinear terms in the van Hove theory 
[19]. The nonlinear term will affect not only the dynamical critical exponent 
but also the power spectrum. The behavior of the power spectrum discussed in 
this paper was almost same as that in the linear response theory. This means 
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Fig. 14. A particle can exchange its energy and momentum by collisions, but it is 
difficult to move far away from the initial position in high density system. 



that the conventional discussions based on the linear response theory may be 
inadequate to discuss the behaviors of the critical fluctuations. Furthermore, 
by using the nonlinear Langevin equation, we can discuss the possibility of 
the dynamical transition in QCD like the ergodic-nonergodic transition that 
is discussed in the glass transition [73] . In the projection operator method, the 
effect of the nonlinearity can be incorporated by introducing nonlinear terms 
as gross variables in the definition of the Mori projection operator. 

In this paper, we have assumed that there is only one gross variable, that 
is, the fluctuations of the order parameter. However, densities of conserved 
quantities also can be gross variables as is discussed in Sec. 3. The most 
promising candidate is the number density of quarks. In the high density 
system, a particle is thickly surrounded by other particles. The energy and 
momentum of a particle is exchanged continuously by collisions, but it is 
difficult to move far away from the initial position, as is schematically shown 
in Fig. 14. Therefore, the number density changes slowly and can be a gross 
variable compared to the energy and momentum densities. This situation is 
probably realized in glass transition where the fluctuations of the number 
density play a role of an order parameter [73]. As a matter of fact, the effect 
of the number density for the critical fluctuations has been studied in the 
hnear response theory [10,11]. 

Definitely, our discussion is applicable to other phase transition like the color 
superconducting phase transition. These subjects are future problems. 
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A Fluctuation- dissipation theorem of second kind 



Sometimes, a Langevin equation is phenomenologically derived, where a noise 
term is artificially introduced to realize thermalization. Then, we assume that 
there exists a relation between a memory function and a noise term like the 
Einstein relation in the Brownian motion. However, in this projection operator 
method, we can prove the exact relation between a noise term and a memory 
term. This is called the fluctuation-dissipation theorem of second kind [38]. In 
the following discussion, we use the Mori projection operator that is defined 
in Sec. 3. 

Let us regard A(x) as a gross variable and derive a Langevin equation for 
A(x). For simphcity, we assume that A(x) does not have a finite value in 
equilibrium, {A{x))eq — 0. Then, the Mori projection is defined by 

P O = (O, A(x')) • (A(x'),^(x"))-^ • ^(x"). (A.l) 
Substituting this expression into the exact TC equation (10), we have 

= (iLyl(x),/lt(x')) . (/l(x'),Xt(x")) ■ /l(x",«) 

t 

+ j ds{iLf{^,s),A\^))-{A{^),A\^'))-'-A{^\t-s) + f{^,t) 

to 

= (iL^(x),At(x')) ■ (A(x'),At(x")) • A(x",t) 
t 

+ J ds(f{^,s),f(^',0)).(A(^'),A''{^"))-'.A(^",t-s) + f(^,t). 

to 

(A.2) 

Here, we put >l(x) = A{x,to), for simphcity. The operator form of the noise 
is given by 

/(x, t) = Qe'^^^'-'°hLA{x), (A.3) 
that has the following properties: 

(/(x,0,A(x')) = 0, (A.4) 
(/(x,i))e, = T:[p/(x,i)] = 0, (A.5) 

where p is given by p = e^'^^/Tr[e~^^]. Here we used that PQ = 0. The 
first correlation means that the noise term is always orthogonal to the gross 
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variable A(x). The second correlation indicates that the thermal average of 
the noise vanishes. 



Now, we introduce the memory function as 



rC 



x,x";t,0) = (/(x,t),/t(x',0)) • (^(x'),^^(x")) 



-1 



(A.6) 



This exact Langevin equation with the Mori projection is, sometimes, called 
the Mori equation. We can see that the memory function is given by the time 
correlation of the noise in the Mori equation. This exact relation is called 
the fluctuation-dissipation theorem of second kind (2nd F-D theorem). In this 
paper, we assume that the exact relation should be satisfied even for the ap- 
proximated memory function. Then, the noise is determined so as to reproduce 
the 2nd F-D theorem, that is, condition (A.6). 



B Derivation of Mori projection operator 

We shall consider the system which is in the thermal equilibrium state with 
small external perturbation {/i^} at t < 0. At t = 0, the external perturbation 
is switched off, and the system relaxes to a new thermal equilibrium state at 
t > 0. The initial density matrix to describe this situation is 



where {A^(x)} is a dynamical variable, whose relaxation we are interested in 
and Zh is the normalization factor. 

The most probable relaxation path at t > in the Heisenberg picture, is given 



At t > 0, the dynamics of the system is governed by the Hamiltonian H. The 
time evolution of the dynamical variable, thus, is defined by 



(B.l) 



by 



{A,{^,t)) = Tr\p{{h,})A,{^,t)]. 



(B.2) 



= e^^*^^(x)e-^ 



(B.3) 
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where ^/i(x) = 74^(x, 0). When the initial deviation from thermal equilib- 
rium is small, the external perturbation {h^} is small and Eq. (B.2) can be 
approximated as 



{A,{K,t))^{A,{K))o 

+ dV/3(A^(x,t), A^i-s!) - (A,(x'))o)o/i.(x'), 

(B.4) 

where (0)o means Tr[p(0)O] and the inner product is the Kubo's canonical 
correlation defined in Eq. (18); 



(F, G)o = y -Tr[p(0)e-^^Fe^^G]. (B.5) 







The small external perturbation {/ijn} in Eq. (B.4) is eliminated using Eq. 
(B.4) at t = 0. Thus, we have 



(A^(x,0> = E / c^VdV'(A^(x,i), A,(x')) • [(A(x'), A(x"))-^],,. • (A^(x 
II II J 



where the inverse is defined by 



E / d=^x[(A(x), A(x'))-^,. • (A.(x'), A(x")) = - x"). (B.7) 

Here, we set (74^(x))o = 0, for simphcity. 

This equation (B.6) can be reexpressed as PA^(x, i) by introducing the Mori 
projection operator. 



PO = E / rfVdV'(0, A,(x')) ■ [(A(x'), A(x"))]-!. ■ A,(x). 

(B.8) 

Namely, the relaxation of an arbitrary operator is extracted by operating the 
Mori projection operator. 
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C The streaming term 



First, it should be noted that the order parameter commutes with the inter- 
action part; 

iLi5a(x.) = 0. (C.l) 
Thus, the streaming term is given by 

e'^^PiLoSa{yi) = e^^* lim P^Sao{x, t) 

= e^" lim J d'^d'^'j^xti^ - x')xo '(x' - x")5a(x"), 

(C.2) 



where 



Xt-t'(x - X') = (5(To(x, t), (5(7o(x', t')) 

r fix 

= J -{ gg(x, t - iX) g g(x^ f ) gg(x, t - iA) | ^(x^ t')} 



u 



EpEp[Ep — Epi 

.2/3+ 

P'^ P P ,(' P _L Vj^, 



+(1 - sin^ - sin^ e;,) ^ ^ 

tjptjpi [tjp + tjp 

(C.3) 



and 



/ ci=^x'xo(x - x')xo '(x' - x") = S^'\^ - x"). (C.4) 

Here, we used the contraction defined in Appendix D. Substituting the above 
expression into Eq. (C.2), it is easy to show that 



lim|xt(x-x') = 0. (C.5) 
Thus, the contribution from the streaming term vanishes. 
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D Contraction of fermion field at finite temperature and density 



In this Appendix, we introduce the contraction of the fermion field. This is 
convenient to implement the calculation with projection operator and we can 
regard it as a correspondence of the propagator in the usual finite temperature 
field theory [45]. 

First, we introduce the Bogoliubov transformation through which the new 
pairs of creation and annihilation operators are defined [74] : 

b^,s ^ cos etBk,, + i sin e+Bl,, (D.l) 
dk,s = cos e^Di,,, + i sin e^Dl^, (D.2) 

where 

sin^ = {ef^'^ + 1)-^ = n^{Ek), (D.3) 

and 



= Vk^ + M2 ±ii = Ek±ii. (D.4) 

Here, M is the mass of the fermion. The above operators satisfy the commu- 
tation relation 



[^k,., Dlj+ = [£'k,., Dlj+ = (D.6) 

All the other commutators vanish. Now, we can define the thermal vacuum 

1^) as 

B^,s\e) = By,,s\e) = D^,s\e) = D^.,s\e) = o. (d.7) 

Then, we can express the statistical average in terms of the vacuum expecta- 
tion value: 

IV[pe,0] = {e\0\9), (D.8) 

where p^q is a density matrix of a thermal equilibrium state. Thus, we can 
use Wick's theorem that simplifies the calculation of the correlation function. 
Wick's theorem tells us that arbitrary correlation functions are expressed as 
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the sum of all possible products of contractions. In our calculation, contrac- 
tions are given by 



go(x^o(x',i') 

go(x',t'>o(x,t) 



g-ip(x-x')gie+(t-t')j._ 



(D.IO) 



E Brownian motion in harmonic potential 



As an exact solvable example, let us consider the one- dimensional classical 
system where a heavy particle with the mass M interacts with heat bath 
composed of many harmonic oscillators [75]. The Hamiltonian of the total 
system is 

The Newton equations of this system are given by 

■ 

Xj = Pj, Pj = -0J%Xj - -^x). (E.3) 

j 



The mass of the particle M is large and the coordinate x and the momentum 
p change slowly in comparison with pi and Xi. In this case, it is natural to 
choose p and x as gross variables of this system. Then, the Mori projection 
operator is defined by 



PG={G,x){x,x) ^x + {G,p){p,p) V 



^(3Mu;',x{G,x) + ^{G,p). (E.4) 

The inner product is given by the canonical correlation for the classical vari- 
ables. 
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{X,Y) = / ^Tr[pe^^Xe-^^r] = Tt[pXY] = (XY), 
J p 



(E.5) 



where p — e '^^/Trfe ^^]. Then, we have the following correlations: 



(p,p) = M//3, 

= = -(p,p)/M, 

{{xi - -fix/uj'^),x) = (p, x) = 0, 
{x,x)^{Mu;',P)-\ 



(E.6) 
(E.7) 

(E.8) 
(E.9) 



When we substitute the Mori projection operator into the TC equation and 
set 0(0) = x{0) and p(0), we have the following Langevin equation 



-x(t) = ^ 



p{t) = -Mujlx{t) - J dsE{t - s)p{s) + fit), 



(E.IO) 
(E.ll) 



where 



1 ^ 7- 

/W = E7acosa;,i(x,(0) - ^a;(0)) + ^p,(0)]. 



a;, 



(E.12) 
(E.13) 



We assume that the frequency of the oscillations distribute continuously 



g{uj) = < 



LO > LOd, 

where uja is a cutoff. Then, the memory function is given by 



(E.14) 



If 7 

^ ^ M J 

" 

sin Udt 

— ^-0 7 — , 

7rt 



(E.15) 
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When wc are interested in the motion with a macroscopic time scale larger 
than 2-11 /ud, the memory function is approximated by 



S(i) ^ 2Eod{t). (E.16) 
Finally, the Langevin equation is give by 



>^f. (E.ir, 

j^p{t) = -Mu'oxit) - Eop{t) + m, (E.18) 
where the noise has the following correlation properties; 



(/W) = 0, (E.19) 
{fit)f{t'))='^E,5{t-t'). (E.20) 



In order to see that the behavior of the memory function strongly depends on 
the choice of the projection operator, we derive the Langevin equation with 
the incomplete Mori projection defined by using one gross variable. 



PG = {G,p){p,p)-'p = ^{G,p). (E.21) 
Then, we have the Langevin equation with one gross variable p, 



d 
dt 

where 



p{t) = -j dsE'it - s)p{s) + fit), (E.22) 



E'{t)=u;l + E{t), (E.23) 
f{t) = -Mu;lx{0) + f{t). (E.24) 

It is clear that E'{t) converges to uJq att ^ oo, although E{t) vanishes. Namely, 
the behavior of the memory function depends on the definition of the Mori 
projection and we must choose the complete set of the gross variables to derive 
the memory function without the long time correlation. 
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